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INVESTIGATION OF TURBULENT FLOW IN A TWO-DIMENSIONAL CHANNEL 1 

By John Laufer 


SUMMARY 

A detailed exploration of the field of mean and fluctuating 
quantities in a two-dimensional turbulent channel flow is 
presented. The measurements were repeated at three Reynolds 
numbers , 12,300, 30,800, and 61,600, based on the half width 
of the channel and the maximum mean velocity. A channel o f 
5-inch width and 12:1 aspect ratio was used for the investigation. 

Mean-speed and axial fl uctuation measurements were made 
well within the laminar sublayer. The semitheoretical predic- 
tions concerning the extent of the laminar sublayer were con- 
firmed. The distribution of the velocity fluctuations in the 
direction of mean flow u' shows that the influence o f the viscosity 
extends farther from the wall than indicated, by the mean velocity 
profile , the region of influence being approximately four times 
as wide. 

Fluctuations perpendicular to the flow in the lateral and 
vertical directions v' and w' , respectively , and the correlation 
u f v* 

coefficient k~ — ~f= were also measured. The turbulent 
■\u' 2 fl v' 2 

shearing stress was computed by three different methods. The 
results show satisfactory agreement for the two lower Reynolds 
numbers. In the case of the highest Reynolds number, however, 
the total shearing stress r obtained from the fluctuation measure- 
ments was approximately 20 percent lower than that computed 
from the mean-velocity and mean pressure measurements. All 
dimensionless mean fluctuating quantities were found to decrease 
with increasing Reynolds number. Measurements of the scales 
of turbulence L y and L z and microscales of turbulence \ u and \ z 
across the channel are presented and their variation with 
Reynolds number is discussed. I sing a new technique, values 
for the microscale \ x were obtained; a new method for estimating 
the scale L x is also given. 

The energy balance in the turbulent flow field wa, s calculated 
from the measured quantities. From this calculation it is 
possible to give a descriptive picture of turbulent-energy diffusion 
in the center portion of the channel cross section. 

For the flow corresponding to a Reynolds number of 30,800 
the energy spectrum of the u' 1 fluctuations at various points 
across the channel, including one in the laminar sublayer , was 
obtained. 

At station y/d=0.4, where y is lateral distance and d is half 
width of the channel, the contribution to the turbulent shear 
stress from various frequency bands was measured and it was 
found that the contribution corresponding to frequencies above 
1500 cycles per second is negligible. Since the spectrum of u' 2 


at this point extends to about 5000 cycles per second it is evident 
that, the high frequencies are nearly isotropic in agreement with 
Kolmogoroff's hypothesis. 

INTRODUCTION 

In recent years a considerable step forward was made in 
the theory of turbulence. Kolmogoroff (reference 1), 
Heisenberg (reference 2), and Onsager (reference 3) ob- 
tained independently an energy-spectrum law that bolds in 
rather restricted types of turbulent flows. This progress 
gave a new impetus to both theoretical and experimental 
investigations. The experimental worker may follow two 
principal methods of approach to the problem. First, he 
may establish flow fields which satisfy sufficiently the 
assumptions of the new theory, namely, that the flow is 
isotropic and of liigh Reynolds number so that the influence 
of viscosity is a minimum and the effect of the turbulence- 
producing mechanism is small. Under these conditions he 
may measure quantities such as correlation functions, scales, 
and microscales that are defined exactly in the flow field 
and may compare his results with those predicted by the 
theory. The main difficulty with this method is to predict 
how closely one has to approximate in the actual flow the 
conditions assumed in the theory. In other words, the 
sensitivity of the theory to deviations from true isotropic 
conditions is not known and in case the measurements do 
not agree with the theory one does not know whether to 
attribute the disagreement to faulty assumptions in the 
theory or to the incomplete isotropic conditions of the flow. 
Furthermore, the imposition of the condition of isotropy to 
fluctuating fields may be a very strong restriction and the 
study of such fields may not yield the complete picture of the 
turbulence mechanism. 

The second method of approach is to establish a simple 
nonisotropic turbulent flow field with well-defined boundary 
conditions and, in the light of the existing theories, to try 
to obtain information on the mechanism of energy transfer 
from the low frequencies of the energy spectrum to the high 
ones. The present investigation of a turbulent channel flow 
has this purpose in mind in its long-range program. 

The difficulties of this method are immediately realized. 
Because of the nonisotropic nature of the flow, characteristic 
quantities such as scale and microscale are no longer well- 
defined and it is very possible that in the first stages of the 
investigation certain quantities will be measured that later 
will prove to be trivial. On the other hand, the main 
advantage of a fully developed channel flow is the fact that, 
in contrast with the flow behind grids, flow conditions are 


i Supersedes NACA TN 2123, “Investigation of Turbulent Flow in a Two-Dimensional 

Channel” by John Laufer, 1950. 


2 


REPORT 1053— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


steady; no decay of mean or fluctuating quantities exists 
in the direction of the flow. Consequently, the turbulent 
energy goes through all of its stages of transformation across 
the channel section — turbulent-energy production from the 
mean flow, energy diffusion, and turbulent and laminar 
dissipation — and one may study these transformations in 
detail. 

It has become clear that the phenomenological theories of 
turbulence, such as the mixing-length theories, have lost 
most of their importance. These theories, developed in the 
late twenties and early thirties, were aimed specifically at 
an evaluation of the mean-velocity distribution in turbulent 
flow. The existing experimental evidence shows clearly that 
the mean-velocity distribution is very insensitive to the 
essential assumptions introduced into the phenomenological 
theories. In fact, purely dimensional arguments generally 
suffice to give the shape of the mean velocity profile with 
sufficient accuracy. For the further development of an 
understanding,of turbulence, detailed measurements of the 
field of fluctuating rather than of mean velocities are neces- 
sary. The program of experimental research of which this 
work is one part is based on this reasoning. It is quite 
apparent and natural that the same conclusions have been 
drawn by workers in the turbulent field elsewhere and that, 
in general, the main emphasis of current experimental inves- 
tigation is the exploration of the field of the fluctuating- 
velocity components. 

The present set of experiments deals with flow in a two- 
dimensional channel, that is, with pressure flow between 
two flat walls. Channel flow of this type is the simplest 
type of turbulent flow near solid boundaries which can be 
produced experimentally. The simpler Couette type flow 
requires that one wall move with constant velocity, a con- 
dition which is difficult to realize experimentally. It can 
be approximated b}- the flow between concentric cylinders, 
but complications due to centrifugal forces arise here. The 
simple geometry of a two-dimensional channel allows an 
integration of the Reynolds equations, and the turbulent 
shearing stress can then be related directly to the shearing 
stress on the surfaces which, in turn, can be determined from 
the mean-pressure gradient or the slope of mean velocity 
profile at the wall. 

The relation between the apparent stresses and the wall 
shearing stress can be used to advantage in two fashions. 
It is here possible to obtain the magnitude of the correlation 
coefficient responsible for the apparent shear by measuring 
only the intensities of the turbulent fluctuations. The 
turbulent shearing stress can also be measured directly by 
means of the hot-wire anemometer. A comparison with the 
shear distribution obtained from mean-pressure-gradient or 
mean-velocity-profile measurements serves then as a very 
useful check on the underlying assumptions, on the one hand, 
and specifically as a check on the reliability and accuracy of 
the direct measurements, on the other. 

Measurements of channel flow have previously been made 
by Doncli (reference 4), Nikuradse (reference 5), Wat ten- 


dor! and Kuethe (references 6 and 7), and Reichardt (refer- 
ences 8 and 9). A few unpublished measurements have been 
made recently at the Polytechnic Institute of Brooklyn . 

Donch’s and Nikuradse’s measurements were concerned 
only with the mean- velocity distribution and are thus of 
not too much interest for a comparison with the present 
set of measurements. Wattendorf measured the intensity 
of the fluctuating- velocity components and then deduced 
the correlation coefficient from the mean-pressure measure- 
ments. The technique for measuring the axial component 
of the velocity fluctuation was well-developed at the time, 
but the cross component was only tentatively measured. 

The most complete set of measurements is due to the work 
of Reichardt, who measured velocity fluctuations in the 
direction of the flow and normal to the wall as well as the 
turbulent shear directly. Reichardt found very good agree- 
ment between the shearing stress determined in these two 
ways; his paper comes closest to the present investigation 
and his results will be used for comparison. The Reynolds 
number in Reichardt’s measurements was 8000, which is 
lower than the range covered by the present measurements. 
A criticism which can be made of Reichardt’s investigation 
is his use of a tunnel of only 1:4 aspect ratio. The two- 
dimensional character of the flow is thus somewhat doubtful. 
Wattendorf’ s experiments were made in a channel of very 
large aspect ratio (18:1) and are thus free of this criticism. 

In a preliminary investigation a channel 1 inch wide and 
60 inches high was chosen. The measurements, however, 
have shown that in this case the scale of turbulence is so 
small that great care must be taken to correct the hot- 
wire readings for the effect of wire length. Measurements 
of the microscales were, in fact, impossible in the 1-inch 
channel. The microscales were about 0.1 centimeter and 
thus smaller than the length of the wire. The corrections 
for the case of measurements of v' and w' were about 30 
percent. Since the method of correction becomes very inac- 
curate for such large ratios of wire length to microscale, the 
measurements were repeated in a 5-inch channel with a 
12:1 aspect ratio. This ratio was still large enough to insure 
two-dimensional flow and the length corrections were greatly 
reduced. 

This investigation was conducted at the Guggenheim 
Aeronautical Laboratory, California Institute of Tech- 
nology, under the sponsorship and with the financial assist- 
ance of the National Advisory Committee for Aeronautics. 

The author wishes to acknowledge the constant advice 
and help in both experimental methods and interpretation 
of results of Dr. H. W. Liepmann during this investigation. 
He also wishes to thank Dr. C. B. Millikan for his contin- 
uous interest in this research. The cooperation of Dr. 
F. E. Marble and Mr. F. K. Chuang is much appreciated. 

SYMBOLS 

x longitudinal coordinate m direction of flow; 

#=0 corresponds to channel exit 
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lateral coordinate; y = 0 corresponds to channel 
wall 

vertical coordinate 
half width of channel (2.5 in.) 
thickness of laminar sublayer 
mean velocity at any point in channel 
maximum value of mean velocity 
instantaneous value of velocity fluctuations in 
direction of mean flow x 
instantaneous values of velocity fluctuations 
normal to mean flow in directions y and z ) 
respectively 

root-mean-square value of velocity fluctuations 
in direction of mean flow x 
root-mean-square values of velocity fluctua- 
tions normal to mean flow in directions y and 
2 , respectively 

pressure at any point in channel 
air density 
total shearing stress 
shearing stress at wall 

friction velocity (V r o/p) 
absolute viscosity of air 
kinematic viscosity 

Reynolds number based on half width of chan- 
nel and maximum mean velocity 
correlation coefficient responsible for apparent 



correlation coefficients as functions of X, Y, and 
Z, respectively 
microscales of turbulence 
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scales of turbulence 

(l z = r Rx dX, L u ”r u d Y, L, = ^R z dz) 

distances (in the x -, y -, and 2 -direction, respec- 
tively) between points at which correlation 
fluctuations are measured 
ratio of square of compensated and uncompen- 
sated fluctuations 

time constant expressing thermal lag of hot- 
wire 

frequency, cycles per second 
fraction of turbulent energy u' 2 associated with 
band width dn 

fraction of turbulent shear u'v' associated with 
band width dn 

dissipation of turbulent energy 
time 

voltage fluctuations 

0 

ANALYTICAL CONSIDERATIONS 


EQUATIONS OF MOTION FOR TWO-DIMENSIONAL CHANNEL FLOW 


The mean- and fluctuating- velocity components are 
denoted by u t and u/ , respectively; Cartesian coordinates 
are designated by .r t ; and p ik denotes the components of 
the stress tensor which includes both the viscous and ap- 
parent (Reynolds) stresses. The Reynolds equation and 
the continuity equation for steady flow can thus be written, 
using Cartesian tensor notation, 
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For channel flow between two parallel flat walls, 
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where x, ?/, and 2 replace x i} and v', and w' replace u/ for 
convenience in writing equations (la) and (lb) in simpler 
form. Equation (lb) is automatically satisfied and equation 
(la) becomes: 


0= ^ ( ~ p - pM ' 2 ' 1 d 


o=A(_,_ p -- ■ > 


)+ 3y( M 3y _pttV ) 

)+ ^( M lr puV ) 


( 2 ) 


In fully developed turbulent flow the variation of mean 
values of the fluctuating quantities with x should be zero; 
that is, the flow pattern is independent of the streamwise 
direction. Hence, equations (2) become: 


1 <*2- 
p dx 


d 2 u du'v' 
' dy 2 dy 


or 


1 dp_ 1 dr 
P dx~ p dy 


1 dp 
p dy~ 


dv /2 

dy 


( 3 a) 


( 3 b) 


Differentiating equation ( 3 b) with respect to x gives 
1 d 2 /? 

p dxdy 

equation ( 3 a) is immediately integrable. Thus 


= 0; consequently dp/dx is independent of y and 


- ^—u'v'+C onstant 

p dx dy 

In the center of the channel the shear vanishes; hence if the 
channel has the width 2d 


and hence 


Constant =- ^ 
p dx 


1 dp du -7—, 

- (y — = u » 

p dx J dy 


Or in nondimensional form: 

y — d dj)__ v du u'v' 

pU 2 dx~~U 2 dy U 0 2 
In terms of the shearing stress at the wall, 
t 0 y — d_ v du u'v ' 


pU 2 


U o 2 dy U 0 2 


( 4 ) 


( 5 ) 


It is evident from equations ( 4 ) and ( 5 ) that t 0 can be 
determined in three different ways: 

(a) From the mean-pressure gradient 


o=-d 


dp 

dx 


(b) From the slope of the mean velocity profile near the 
wall 

du > 


T 0— P l -7— 1 

\dyJv = 0 

(c) From a direct measurement of u'v' 

( -7—7 du\ d 

To= ( puv 


If y/d is not too small, 


du , , -7-7 

IXJ - << P UV 


and then 


t 0 = pu'v' 


y—d 


The technically simplest way to determine r 0 , and thus in 
fact the complete shear distribution, is a measurement of 
dp/dx. This is the method applied in most investigations. 

To determine the slope of the velocity profile near the wall, 
the profile has to be known up to points very close to the 
wall, that is, at least to distances y/d= 10“ 2 . This in general 
requires the use of the hot-wire anemometer and very 
precise measurements. 

The third method requires a direct measurement of the 
correlation between the axial and lateral velocity fluctuations. 
The technique of this type of measurement is known .and 
was first applied by Reichardt and by Skramstad and, in 
somewhat different form, in recent investigations at the 
National Bureau of Standards, Polytechnic Institute, of 
Brooklyn, and California Institute of Technology. 

In the present investigation all three methods have been 
applied and the results compared, with the exception of the 
flow at the lowest Reynolds number. In this case the 
pressure gradient is extremely small (approx. 0 .C 003 mm of 
alcohol/cm) and reasonably accurate measurements were not 
possible. This comparison of the three methods has the 
advantage that it gives a good indication of the absolute 
accuracy of measurements of the fluctuating quantities and 
the correlation coefficient k. 

ENERGY EQUATION 

Writing equation (la) in the form 


dUi dp . 


(6) 


and multiplying it by u t , the following relation is obtained: 

1 dUjUjU k _ dpUj , d 2 Uj 

2 p dx k dXj ' pUl dxjdxj 

Transforming the last term by partial integration the energy 
equation of the mean flow becomes 
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1 dUjUjUk djnij 1 d 2 UjUj /duA /duA . 

2 p dx, : dxi 2 p dxjdxj M \d^/ \ dx / 

The velocity perturbations can now be introduced: 

u x =u-{-u' 
u 2 = v ' 
u 3 =w' 


Since the velocities are independent of the coordinates x and z 
the following equation is obtained after averaging: 


—r~,du . 
pu v j — r pu 
ay 


du'v' . 1 (I 
~dy ^ 2 p dy 


v '(u' 2 + v ,2 + w ,2 )=-u 

dx 


<>v'V ■ 

dy ^ 


see whether experimental results confirm this relation. Taking 
as an example results from the measurements at I\= 30,800 
and y/d = 0.5 (where q 3 /D ~ 11 X 1C 3 and vq 2 /X 2 ^ 1 .7X 1 0 3 ) 

q = 52 centimeters per second 

\ v =0.5 centimeter 

therefore 

qX 2 52X0.25 
vD 0.155X12.7 

This ratio seems to be fairly constant for different Reynolds 
numbers and across the channel cross section. 

In view of the above dimensional considerations, the 
energy equation may be written 


d 2 u 1 

1 


d?_ 

dy 2 


(u' 2 + v' 2 -\-w' 2 )— n 


/du/\ /du/\ 

\dxjJ\dxj/ 


du d ( ,u'~-\-v' 2 +w ' 2 

T ( brddiV v 2 



(8) 


Making use of equation (3a) this simplifies to 


dy dy 


d 2 


2 M dy 2 ' " 


The equation expresses the fact that the energy produced 
by the turbulent shear forces at a certain point is partly 
diffused and partly dissipated. This equation has been 
used in estimating the energy diffusion across the channel, 
since from the measured quantities the production term and 
dissipation can be calculated. 


This form was obtained by Von Karman (reference 1C) while 
discussing a nonisotropic flow in terms of the statistical 
theory of turbulence. In order to see the relative orders 
of magnitudes of the different terms in the equation, Von 
Karman expresses them in terms of a single velocity 
q=-yj l u r2 -\-v /2 -\-w r2 and a characteristic length I) correspond- 
ing to the width of the channel. Since 


-«-«V=0(g 2 ) 

p 



the above equation may be written as 


,1, f rAi £+A 3 $-+A< 


"2 2 , j v jl 
D 2 ^ 5 X 2 


where the coefficients have the characteristics of correlation 
functions and X is the microscale of turbulence. For a high 
Reynolds number flow qD/v^> 1, it follows that vq 2 jD 2 <^(flI). 
Thus the second and third terms on the right side of the 
equation may be neglected. Since the pertinent quantities 
have been measured during the present work, the order of 
magnitude of these terms can be directly evaluated and the 
omission of the terms is found to be justified. The above 
equation thus contains only terms of the forms q z /D and 
vq 2 /\ 2 ; these terms should be of the same order of magnitude 
and therefore q\ 2 /vD= 0(1). It is of considerable interest to 


EQUIPMENT ANI) PROCEDURE 

WIND TUNNEL 

The investigation was carried out in the wind tunnel 
shown in figure 1. The turbulence level is controlled by a 
honeycomb and seamless precision screens, followed by an 
approximately 29:1 contraction. The screens have 18 
meshes per inch and a wire diameter of 0.018 inch. The 
honeycomb consists of paper mailing tubes, 6 inches long 
and 1 inch in diameter. 


Wind direction .. -Screens (! 8 mesh) 



5" channel ^ 

Access doors-*'' ''Honeycomb 

(Tubes 6" long and /" diam.) 


A 

□5a 


Wind direction ^ Box section \ 



Counter vanes-- 
Conical diffuser S j 


|j=£ 


Propel tens-*'' 


Figure 1 .— Diagram of two-dimensional tunnel. 


The over-all length of the channel is 23 feet. At the en- 
trance section it is 3 inches wide and has an aspect ratio of 
20: 1 . In a distance of 7 feet it expands to a width of 5 inches 
and the aspect ratio is reduced to 12:1 . The walls of the exit 
portion of the channel (about 6 ft) are made of %-inch- 
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thick plywood with a %-inch birch inside cover. They were 
specially treated in order to acquire a smooth finish and were 
reinforced in order to avoid warping (fig. 2) ; in spite of this a 
few percent of width variation existed 

NOT RFPRODUCIBLE 


Figure 2. — Exit of channel. 


The tunnel is operated by a 62-liorsepower stationary 
natural-gas engine — normally operating at a fraction of its 
rating — which drives two eight-bladed fans. The speed is 
remotely controlled by means of a small electric motor which 
drives the throttle through a gear and lead-screw system. 

The experiments were carried out at speeds of 3, 7.5, and 
15 meters per second. 

TRAVERSING MECHANISM 

Figure 3 shows the type of traversing mechanism used 
during the experiments. It consists simply of a micrometer 
screw on which the hot-wire support is fastened. The sup- 
port can rotate freely in a plane perpendicular to the air flow; 
thus, the hot-wire may be adjusted exactly parallel to the 
wall. 


/ 

// 

- 


Figure 3.— Traversing mechanism. 

The zero reading of the traversing mechanism (y = 0) was 
carefully found using the following method: The hot-wire was 
placed close to the wall (approx. 0.025 cm away) ; the distance 
between the wire and its image in the polished wall was 
measured by means of an ocular micrometer. The position 
of the wire is, of course, one-half of the observed distance and 
could be determined with an accuracy of ±0.0005 centi- 
meter. 

In order to obtain the pressure distribution along the 
middle of the channel a 6-inch-long thin-walled tube %-ineh 
in diameter was used with a small static-pressure hole drilled 
close to its end. The tube was free to slide in a support at 
the entrance of the channel so that the position of the pressure 
hole relative to the channel exit could be changed. The tube 
was kept straight and under tension by means of weights. 

HOT-WIRE EQUIPMENT » 

All velocity measurements were made with hot-wire ane- 
mometers. The frequency response of the amplifier-compen- 
sator unit of the hot-wire equipment, using a 0.00024-inch 
wire at standard operating conditions, is fiat from approxi- 
mately 2 to 10,000 cycles per second. 

The compensation of the hot-wire for thermal lag is ac- 
complished by a capacity network. The range of time 
constants was chosen from 0 to 1 millisecond corresponding 
to the characteristics of 0.00005- to 0.00024-incli wires at the 
operating conditions employed in general at GALCIT. No 
attempt was made to extend the range of compensation to 
larger values of M, since in this case the noise level soon 
becomes appreciable. 

The correct setting of the compensating unit was found 
using the square-wave method described by Kovasznay (ref- 
erence 11). The time constants of the wires used for turbu- 
lence measurements fell between 0.1 and 0.9 millisecond. 

The output readings were taken with a thermocouple and 
precision potentiometer. 

Mean-speed measurements. A 1/2-mil (0.0005-in.) plat- 
inum wire of 1 -centimeter length was used to measure the 
mean speed. The measurements were made by the constant- 
resistance method. This method has the advantage of 
keeping the wire temperature constant through the velocity 
field. 
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Turbulence measurements. For the investigation of tur- 
bulent fluctuations, 0. 00024-inch Wollaston wire was used. 
The wire was soft-soldered to the tips of fine sewing needles 
after the silver coating had been etched off. The two lateral 
components v' arid w' and the correlation coefficient 

k= j - — = were measured using the X-wire technique. The 

VW *’' 2 

method was essentially the same as described in detail in 
previous reports from this laboratory. (See, for example, 
reference 12.) ■ c . 

The X -meters for shear measurejnents had angles between 
the wires of approximately 90°. The angles of the v',w'~ 
meter were of the order of 30°. The wire length of the u'- 
meter was about 1.5 millimeters and that of the v r -meter 
was 3 millimeters. 

The parallel-wire technique used recently at the Polytech- 
nic Institute of Brooklyn was also tried in order to obtain 
the lateral component of the velocity fluctuation and the 
shear close to the wall. The parallel-wire instrument is 
superior to the X -meter since it allows an exploration of the 
turbulent field up to distances of a few thousandths of an 
inch from the wall. However, it was found impossible to 
obtain reliable values with this instrument. The correc- 
tions due to w' -fluctuations and to unequal heating of the 
wires were pronounced and not easily accountable. The 
method was therefore temporarily abandoned after consid- 
erable time and labor had been spent on it. 

Measurements of correlations between two points. The 
correlation functions between values of u' at points along 
the y- and 0 -axis, that is, 

R _ u'{0)u'(Y) 

^u\oy^u f (jY 

v/(0 )u'(Z) 

Ji 2 — — = 

z / o 


were measured using the standard technique (reference 12). 
The scales of turbulence L u and L 2 and the microscales \ u 
and X 2 at different points across the channel were obtained 
from these measurements. 

Measurements of X x . — A method suggested by Townsend 
(reference 13) was applied to measure X x . Using an electronic 
differentiation circuit the amplifier output signal was dif- 
ferentiated and X x computed from the following relation: 



of the ^'-fluctuations. From these counts X x may be cal- 
culated directly by assuming a normal and independent 
distribution for both u' and du'/dx: 


1 7 r 

r-=-jyX Average number of zeros of u' per second 
Xx ' 0 


It is known that the distribution of u' is closely a Gaussian 
one even in nonisotropic turbulence (see, for instance, refer- 
ence 15); however, for the case of du'/dx a small deviation 
from the normal distribution was found (reference 13). For 
the preliminary measurements of X x reported presently, no 
corrections were applied as yet for this effect. Figure 4 
shows an oscillograph trace of the ^'-fluctuation in the mid- 
dle of the channel at 7^ = 30,800. The trace represents an 
interval of approximately 1/20 second. 



Figure 4.— Typical oscillogram of u '-fluctuation at /?=30,80() and y/d=l.Q. Horizontal line 
corresponds to u'(t)= 0; interval is approximately 1/20 second. 


Measurement of L x . — The following simple procedure 1 
was applied to obtain a rough estimate of scale of turbulence 
corresponding to correlations between points along the x - 
axis: Denote by F~ 2 ( 11 ) the fraction of turbulent intensity 
which is contributed by frequencies between n and n-\-dn ; 
that is, 

u f (n)- dn=u ,2 F~' 2 ( 71 ) dn 

and thus 



F lL '*(ri) dn= 1 


Consider now an uncompensated hot-wire. If the time 
constant of this uncompensated wire is M, the response 
will be 

\u (jiy\uncomp. 

where 


[tt (ft y\uncomp. 


u'(n) 2 

T+M 2 n 2 


The error involved in the approximation yu — can be 

estimated and it may be seen that the above relation holds 
with reasonable accuracy over the large center portion of 
the channel. 

A new r technique for the measurement of X x has also been 
applied. The method is described in detail in reference 14 
and consists in counting the zeros of an oscillograph trace 


The total intensity for the uncompensated wire will then be 
given by 3 


(u rl ) 


tin comp. 



FTjKn ) 

1 + M 2 ri 2 


dn 


2 This method was suggested by Dr. H. W. Liepmann. 

3 Formulas of this general type have been proposed by Kamp6 de Ferict and by Frenkie 1 
for determining the spectrum of turbulence from uncompensated hot-wire measurements by 
varying M (reference 16). 


981431 — 52 2 
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Thus, if the ratio of the response of the same wire compen- 
sated and uncompensated is denoted by < 7 , 


1 c m Fu*(n) J 

~J. T+MV in 


(9) 


In order to estimate L x , F u f 2 (n) is assumed to have the simple 
form 

x 2 L X \U 0 
Fu<n ) — _ j 2 


1 +tt 2 


L 2 

Uo 2 


which was given by Dryden (reference 17) and corresponds 
to an exponential correlation curve. 

Then 



or with 

V == LxH jUg 


then 


Hence 


a= MUJL X 



drj 1 

(\+T)*)(l-\-a 2 T) 2 ) 1-fc* 


L x 


MU g 

a — 1 


Thus by measuring the ratio of the mean squares of the 
fluctuating velocities with and without compensation, L x can 
be estimated. 

The procedure was checked in the flow behind a Jo-inch 
grid where the lateral scale was measured and where because 
of isotropy the relation L x = 2 L y should hold fairly well. 

The results were: 

<r= 1.85 

M=7X 10~ 4 second 
l ' 0 = L5X 10 3 centimeters per second 


therefore 

L x = 1.24 centimeters 

Direct measurements gave L v = 0.553 centimeter. Hence 
L x =l.ll centimeters which shows that this method of 
estimating the value of L x is satisfactory. 

Measurement of turbulent-energy spectrum. — In order to 
be able to measure accurately the energy distribution over a 
wide frequency band, it was found necessary to use extremely 
thin wires with small time constants and thus to increase the 
signal-to-noise ratio. Wires of 0.00005-inch diameter, op- 
erating with time constants of approximately 2 X 10 -4 
second, were used. In this way it was possible to detect 


energy values in the high-frequency bands that were 10~ 7 
times the value corresponding to zero frequency. 

The hot-wire signal was fed into a Hewlett-Packard wave 
analyzer, the output circuit of which was somewhat altered 
in order to be able to feed the output signal into a root- 
mean-square voltmeter. The averaging characteristics of 
the root-mean-square meter were checked and, for a narrow- 
band-width signal, were found to give the same mean values 
as those calculated from thermocouple readings. 

With this method the energy of the i/'Mluct nation was 
measured at different points across the channel. At y/d= 0 A 
the spectrum of the turbulent shear u'v' was also obtained. 4 
An X-type meter was used, both wires making an angle of 
45° with the free-stream direction. If e 2 and e 2 are the mean- 
square voltage fluctuations (corrected for time lag) across 
the two wires it may be shown (reference 12) that 

ei 2 — e 2 = Cons tan t X u' v' 

If 6 ?i and e 2 are fed into the wave analyzer then the above 
relation holds for all band widths dn; that is, 

e i (r) 2 — e 2(71 ) 2 = Const an t X u' v\n) 

PRELIMINARY MEASUREMENTS IN A 1-INCH CHANNEL 

In the initial stages of the turbulent-channel-flow investi- 
gation a two-dimensional channel of 1 -inch width was used. 
Measurements of mean and fluctuating velocities and of the 
correlation coefficient k were completed. The scales L y and 
Lg were also measured at the channel center and were found 
to be about 0.2 to 0.3 centimeter. This small-scale turbulence 
existing in the channel imposed a definite limitation on .the 
accuracy of the fluctuating measurements. Wire-length 
corrections as high as 30 percent had to be applied to the 
measurements of v' and w'. Furthermore, no microscales , 
could be measured accurately; thus one of the objectives of 
the investigation, the calculation of the energy dissipation ! 
across the channel, could not be obtained. Nevertheless, 
results show good consistency; the three independent meas- , 
urements of the turbulent shear indicate satisfactory , 
agreement. I 

It is of interest to present these preliminary measurements, , 
the Reynolds number of which was 12 , 200 , and to compare j 
them with those obtained in the 5-inch channel. Figure 5 | 

shows the distribution of all the measured quantities in the 
1 -inch channel and it may be directly compared with figure 10 
which shows the quantities measured in the 5-inch channel at f 
nearly the same Reynolds number. 

Although the slopes of the mean velocities at the wall are 
almost the same in each case (for the 1 -in. channel r 0 /pL r 0 2 = 1 

0.0019; for the 5-in. channel rjpl Y= 0.001 8), the velocity | 
ratios ujV 0 in the 5-inch channel seem to be lower across 
most of the channel. This indicates that a more intense j 

* This method appears to have been first used by I)r. Stanley Corrsin in a turbulent jet j 
(private communication). 


INVESTIGATION OF TURBULENT FLOW IN ^ TWO-DIMENSIONAL CHANNEL 



x/0 4 


x/0 4 


^urbulent-enorgy production takes place in the 5 -inch 

channel. Indeed, the turbulent- velocity fluctuations u'/( 
v'jl o, and w'/l'o (especially u'/U 0 ) increase at a faster rate 
toward the wall in the wider channel, although their values 
at the channel center check very well in the two cases. The 
shear distribution r/pf// is almost the same for the two 
channels; it follows, therefore, from the relation 

uV , v/_ v' 

Uj Uo Uo 

that k must have a higher absolute value in the 1-inch 
channel since both u'/l \ and v' / 1 \ are smaller there. Indeed 
the maximum value of k in this case is ~0.6, while in the 
5 inch channel it is — 0.5. 

As a matter of interest it could be mentioned that if the 
basis of comparison is not Reynolds number but maximum 
mean speed (for the 1-in. channel / 0 =15 in/ sec) then figure 
5 shows a distribution of u' , r', and w' very similar to that 
of figure 18, the absolute values being somewhat lower in 
the 5-inch channel. 

RESULTS AND DISCUSSION 

MEAN- VELOCITY DISTRIBUTION 

A careful study of the two-dimensional nature of the 
channel flow was first made. Mean- velocity measurements 
were carried out at approximately x=—2 inches at different 


heights in the channel: At positions 6 inches from the 

bottom, 6 inches from the top, and at the middle. Agree- 
ment among the three sets of measurements confirmed the 
two-dimensionality. A further check was made on possible 
end effects that might influence the results. The length of 
the channel was extended by 6 inches and the mean-velocity 
measurements were repeated at x=— 8 inches. No change 
in the profile was noticed. Figure 6 shows the mean-velocity 
distributions at three Reynolds numbers, 7?=61,600, R= 
30, 800, and R= 12,300. The distributions follow Von 
Karimin’s logarithmic law very well, except, of course, near 
the wall and at the center of the channel (fig. 7). The 
values of l \ were obtained from the velocity gradients at 
the wall. 



Figure 6.— Mean-velocity distribution. Measured points close to wall are not shown. 
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Measurements were made with special care close to the 
wall. Velocities at a large number of points were recorded 
within the laminar sublayer in order to establish with 
reasonable accuracy the shape of the velocity profile at the 
wall (fig. 8). The thickness of the laminar sublayer (the 
point where the velocity distribution deviates from the 
logarithmic law) was found to be 5«30 r/f r (fig. 7). 



In figure 8 it can be seen that for the case of the lowest 
Reynolds number a few points near the wall indicate very 
low velocities. Since the point at which yld=0 has been 
determined with great accuracy and since the equations of 
motion of the channel flow require a negative curvature for 
the velocity profile, it was concluded that the hot-wire 
indicated too low velocities near the wall for the lowest 
Reynolds number. (The dashed curve for this case in fig. 8 
shows the interpolated velocity distribution for0<^t//d<\0.01 .) 
It can be shown that very high local- velocity fluctuations 
cause the hot-wire to read velocities lower than the true 


mean value. Since the mean current of the wire varies with 
the fourth root of the velocity, 


Vu=Y 


T/ , w'Y , /VY 

! , /VYT /8 

LO-hiO+U) 

+U)J 


where u is the velocity to which the hot-wire responds and 

u is the true mean velocity. Near the wall ^4-^ 

and can therefore be neglected. Because of the cosine-type 
directional characteristics of the hot-wire the effect of the 
w'-fluctuation can also be neglected. The velocity measured 
by the wire can be written, after expanding the above 
expression, as 


u = u 


[>-£(*)■+-] 


Since the average of the odd terms is zero the series con- 
verges rapidly. In the region in question, the velocity 
fluctuations obtained were very high indeed (S'/^^>0.30). 
Their absolute values, however, are probably even larger, 
since the nonlinearity effect of the large fluctuations on 
the hot-wire response still further amplified by the very low 
mean velocities is not taken into account because the wire- 
length corrections are neglected. The mean- velocity cor- 
rection therefore given by the above relation is too small to 
give the required negative curvature of the profile. 

TURBULENCE LEVELS 

Figures 9 to 14 represent the results of measurements of 
the three components u', v', and w' of the turbulent-velocity- 
fiuctuation distribution in the channel for the three Reynolds 
numbers. 



t 






I 
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Figure 11. — Velocity fluctuations u'. 





Figure 14.— Distribution of velocity fluctuations w'. 

The velocity fluctuations u ' relative to local speeds 
increase very rapidly near the wall as is shown in figures 9 
and 10. Measurements very close to the wall indicate that 
u'/u reaches a maximum within the laminar boundary layer 
(yUr/v^ll) and it tends toward a constant value at the wall 
which is independent of the Reynolds number. This point 
will be discussed in detail under “Reynolds Number Effect. ” 

The absolute values of the distribution of u ' show the same 
general shape as that obtained by Reichardt (reference 9), 
having the characteristic maximum near the wall and thus 
showing the strong action of viscosity even for values of 
2 / — 45 . 

Using the X-type hot-wire technique for obtaining the 
velocity-fluctuation components v' and w', no measurements 
could be obtained near the wall. Figures 13 and 14 show 
that, while in the center of the channel the magnitudes of 
v f and w' are the same, w' increases faster toward the wall. 
This agrees with the ultramicroscope measurements in a pipe 
by Fage and Townend (reference 18). 

No length corrections were necessary to the measurements 
of u' except near the wall; however, no corrections were 
applied in this region since no measurements of \ z could be 
made. In this region, furthermore, the fluctuations are 
very large and the values given in figure 10 must be accepted 
with reserve. The hot-wire response for large velocity fluc- 
tuations is not well- understood yet and no correction was 
attempted. Length corrections were applied to the measure- 
ments of v ' and w ' . 

CORRELATION COEFFICIENT ANI) SHEAR DISTRIBUTION 

The correlation coefficient is fairly constant across most 
of the channel (fig. 15) as indicated already by Wattendorf 
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Figure 15.— Correlation-coefficient distribution. 


and Kuethe (reference 6). The maximum values of k ob- 
tained decrease slowly with increasing Reynolds number and 
thus show a definite Reynolds number dependence. Existing 
results indicate, however, that parameters other than Reyn- 
olds numbers also influence the value of k. The following 
table shows magnitudes of the maximum correlation coeffi- 
cient as obtained by different investigators working with 
channels of various widths and with various Reynolds 
numbers: 

Channel width 


Experiments by — 

(cm) 

R 

k m ax 

Reichardt (reference 9) 

24. 6 

8, 000 

-0. 45 

Laufer (present paper) 

2. 5 (1 in.) 

12, 200 

-. 63 

Laufer (present paper) . . 

12. 7 (5 in.) 

12, 300 

-. 50 

Wattendorf (reference 7) 

5. 0 

15, 500 

-. 52 

Laufer (present paper) 

12. 7 (5 in.) 

30, 800 

— . 45 

Laufer (present paper) __ _ _ 

12. 7 (5 in.) 

61, 600 

-. 40 


Wattendorf obtained his value of k from his measurements 
of the quantities u'/U 0) v'/U 0 , and u'v'IU 0 2 . Unfortunately 
his preliminary measurements of v f are incorrect; their 
magnitude is larger than that of u' at all values of y across 
the channel contrary to results obtained by Reichardt and 
by the present writer. Wattendorf himself points out the 
inaccuracy of the values of v' in his paper. Because of the 
large values of V his computed correlation coefficients are 
too low. The value k max =- 0.52 listed in the above table 
was obtained by using Wattendorf^ values u'r'IU 0 2 = 0.001 
and u/U 0 = 0.055 and the value v f /U o =0.035 (y/d= 0.5) 
obtained by the author for both the 1-inch (2.5-cm) and 
5-inch (12.7-cm) channel at Ii= 12,200 and R= 12,300, 
respectively. 

The variation of k indicates that for the same Reynolds 
number the correlation coefficient tends to increase with 
increasing maximum mean velocity (i. e., in a channel of 
decreasing width). 

Figures 16 to 18 show the distribution of all the measured 
quantities including the shear distribution. For the case of 
the two higher Reynolds numbers, three independent 
measurements were made for the determination of the shear 
distribution by methods indicated in “Equations of Motion 
for Two-Dimensional Channel Flow.” Figure 19 indicates 
that consistent results are obtained for the value of r whether 
calculated from the pressure gradient along the channel or 
from the velocity gradient at the wall. The turbulent-shear 
distributions obtained from hot-wire measurements and 
those calculated from figure 17 show satisfactory agreement. 
However, for the flow at the highest Reynolds number, the 


hot-wire measurements gave a 20-percent-lower value for 
the shear coefficient, as may be seen in figure 18. 5 At 
present, the writer can give no satisfactory explanation for 
this discrepancy. In the flow at the lowest Reynolds 
number, no pressure measurements were made because the 
very low pressure gradient (approx. 0.0008 mm of alcohol/cm) 
did not permit accurate measurements. However, in this 
case the wall shear computed from the measured mean 
fluctuating quantities u' } v', and k and the mean-velocity 
measurements may be compared with that obtained in the 
1 -inch-channel flow having nearly the same Reynolds 
number. The comparison shows satisfactory agreement: 
For the 1-inch channel T o /pU o 2 =0.00\9 ; for the 5-inch 
channel T u /pU o 2 =0.00\8. 

SCALE ANI) MICROSCALE MEASUREMENTS 

For a further understanding of the structure of the tur- 
bulent field, correlations of ^'-fluctuations at two different 
points were carried out. Since the field is not isotropic, 
the scale and microscale measurements were repeated both 
in the y - and ^-direction for different values of y/d. 

/^-correlation. — Figure 20 shows typical /^-correlation 
distributions at different values of y/d corresponding to 
//= 30,800. For larger values of z , inaccuracies in the 
measurements did not permit the exploration of the negative 
region of the correlation distribution. From the measured 
distributions of R z at four stations across the channel and 
for different values of /?, the values of L z and \ z were cal- 
culated (figs. 21 and 22). In these figures L z is seen to be 
decreasing uniformly toward the wall, while \ z reaches a 
definite maximum at about y/d= 0.7 and then decreases 
with decreasing y/d. 

R v correlation. — Distributions of the /^-correlation at 
a given value of y/d were obtained by fixing the stationary 
hot-wire at the given value of y/d and traversing with the 
moving hot-wire away from the fixed one toward the channel 
center. Values of L u and \ v were then calculated from these 
distributions of R v . In the region 1 .O^y/d^O.l the grad- 
ients of the various quantities are slight; no significant 
asymmetry in the function R v is therefore expected. Meas- 
urements of the top part of the /^-correlation distribution 
(1 R v ^>0 .8) by Prandtl and Reichardt (reference 19) 
justify this belief fairly well. 

By comparing figure 20 with figure 23, the strong effect 
of the existing nonisotropy can immediately be seen. Al- 
though both R v and R z correspond to Von Karman’s g - 
function (reference 10) they exhibit different behaviors 
across the channel. The R z - correlation function falls more 
and more rapidly to zero with decreasing y/d. Traversing 
from y/d — 1.0 to y/d— 0.70 the distribution of R y is found to 
behave similarly; however, the dashed curve measured at 
y/d=0A (fig. 23) indicates considerably increased correla- 
tions for larger values of Y. It is thus seen that around 
y/d= 0.7 there is a definite decrease in energy content of the 
fluctuations having low frequencies. Further consideration 
of this fact is given later in the discussion of the energy 
balance in the channel. 


5 This result was accepted after the absolute values were carefully checked; the />'-wire 
response was compared with '-fluctuations in an isotropic field and the two-dimensional 
character of the flow was checked. 


I 


INVESTIGATION OF TURBULENT FLOW IN A TWO-DIMENSION AL CHANNEL 


13 



V/d 

Figure 16. — Distribution of mean and fluctuating quantities. 72=12,300. 
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Figure 17.— Distribution of mean and fluctuating quantities. 72=30,800. 
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Figure 18.— Distribution of mean and fluctuating quantities. #=61,600. 
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(a) #=61,600. 

(b) #=30,800. 

Figure 19.— Pressure distribution along channel. Solid lines indicate pressure gradient 
obtained from measured du/dy. 
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Figure 20— #.-correlation. #=30,800. At y/d= 1.0, £, = 1.65 centimeters, X,=0.70 centi- 
meter; at y/d= 0.7, £, = 1.55 centimeters, X,=0.81 centimeter. 
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The values of L y = j R y dY (fig. 21) are not so reliable as 

those of L z since R v could not be measured at sufficiently 
large values of Y because of the limitation of the traversing 
mechanism used. 

The distribution of \ y across the channel is similar to that 
of X 2 (fig. 22). From y/d= 1.0 to y/d~ 0.7 both \ y and X 2 in- 
crease almost proportionally to u' , indicating that the tur- 
bulent-energy dissipation W oc u /2 /\ 2 is approximately con- 
stant in this region. It should be noticed, however, that 
X 2 is consistently larger than \ y throughout the channel cross 
section. Some measurements of X /; close to the wall are also 
indicated in figure 23. No length correction was found to 
be necessary for the measured values of \ y and X 2 . 

A rough estimate of L x by the method already described 
gave a value approximately twice that of L z at the center of 
the channel (fig. 21). Its value increases to a maximum at 
y/d~ 0.5 and then decreases rapidly. No values for L x are 
given for the lowest Reynolds number; in this case the value 
of a is very close to unity and, since L x is proportional to 

— - — , a should be known within an accuracy of 1 percent to 
a — 1 

give consistent results for L x . Unfortunately measurements 
of L x cannot be made with this accuracy, particularly when 
the ^'-fluctuations are of rather low frequencies as is the case 
for 7^= 1 2,300. It should be mentioned that the accuracy 
of the determination of L x is more limited by the inaccurately 
measured value of a than by the fact that an approximate 
spectrum function F~z(n) is used in equation (9). 
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clue to the large noise- to-signal ratios. 


As a matter of interest the results of some measurements of 
\ x obtained with the zero-counting method are compared 
with those described above. (See fig. 22.) Since the validity 
of the assumptions involved in this method is not clarified 
yet, these measurements should be regarded as preliminary. 

SPECTRUM MEASUREMENTS 

The spectra of the velocity fluctuations u' 2 have been ob- 
tained at various values of y/d across the channel. With 


improved instrumentation it was possible to reduce experi- 
mental scatter by a considerable amount. The accuracy of 
the present measurements is believed to be within ±10 
percent with the exception of values corresponding to low 
frequencies (n<100 cps) because of the large-amplitude 
fluctuations and with the exception of values corresponding 
to high frequencies (#>4000 cps) because of noise and of 
possible wire-length effect. A typical spectrum distribution 
taken at y/d= 1.00 is shown in figure 24 . (The measured 
spectra at various positions of y/d are given in table I.) No 
attempt is made to compare the measured distributions with 
existing theories since the restrictive assumptions of these 
theories (isotropy and flows at high Reynolds numbers) are 
not satisfied in the present experiments. The only frequency 
range where comparison is possibly justified is the viscous 
region, that is, the high-frequency part of the spectrum where 
viscosity plays a dominant role. Unfortunately here the 
accuracy of the measurements is not good enough to afford 
any definite conclusions. Therefore, the rather close agree- 
ment of the measured spectrum (for all values of y/d except 
in the laminar sublayer) with the /r 7 -law predicted by Heisen- 
berg (reference 2 ) should be accepted with reserve (fig. 24 ). 




Figure 23.— # y -correlation. .R=30,800. At y/d= 1.0, L y = 1.3 centimeters, X y =0.50 centi- 
meter; at y/d= 0.7, L y = 1.2 centimeters, X v =0.56 centimeter. 




J 
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Fin i re 24. — Frequency spectrum of u'- at channel center. Scale values within parentheses 
correspond to curve on left. 



For a comparison of the energy spectra at different points 
in the channel, figure 25 shows the distribution of )i 2 F U ' 2 (n) 
instead of that of F^(n) as a function of t lie frequency. 
This method of presentation emphasizes energies correspond- 
ing to the higher frequencies but gives a clearer over-all 
comparison of the different sets of measurements. The 
characteristic maximum of the distribution of \ x immediately 
becomes apparent in this figure. From the fact that the 


TABLE I.— MEASURED VALUES OE THE SPECTRUM 


OF u ' 2 





F~(n) 



Fre- 

quency 

(cps) 



(sec) 



1.0 

5 = 0.7 

§=0.4 

§=0.1 

=0.005 


d 

d 

d 

d 

d 

30 






(772 X 10-5 
1786 





I 

40 

561X10- 5 

585 X lO- 5 

571 X lO" 5 

664 X 10-5 

\ 

(618 

1604 

50 





519 

60 

349 

464 




436 

(373 

1345 

70 



\ 





1 

80 

273 

211 

252 

308 

309 

90 





/254 

1272 





100 

125 

186 

141 

182.5 

226 

\ 

222 

(157 

11.50 





1 

150 

175 

95.4 

77. 3 

91.6 

140 

1 

j 

(109 
1114 
( 86. 5 
( 74.9 



| 



200 

68.2 

52.6 

56. 8 

97.2 


( 63.6 
1 63. 6 

250 

23. 1 

37.9 

45.8 

64.9 

1 

( 38. 6 
1 38.0 

300 

( 33.9 

23.5 

32.2 

44.8 

23.8 

l 34.6 

25.6 

34.7 

51.3 

21.5 

400 

19.8 

15. 0 

19. 17 

25.2 

8.03 

500 

18.9 

8. 09 

12. 25 

15.5 

3. 74 

600 

7.57 

5. 86 

7. 6 

12.5 

1.83 

700 

5. 27 

4.0 

5. 14 

8. 73 

1.06 

800 

3. 31 

2. 67 

3. 53 

6.07 

.531 

IKK) 

2.29 

1.77 

2.51 

4.04 

.352 

1000 

1.50 

1.0S 

1.76 

2. 99 

.185 

1200 

.764 

.677 

1.01 

1.80 

.090 

1500 

.295 

. 257 

.379 

.787 

.016 

17.50 

. 147 

. 16 
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correlation coefficient is a Fourier transform of the spectrum 
function it follows that 



1 _4tt 2 r 
\ x Jo 


n 2 F u ,2 ( n )dn 


Thus \ x can be computed from the distributions of 
n 2 Fu'*(n). The calculated values check within a few per- 
cent with direct measurements (fig. 22). 

REYNOLDS NUMBER EFFECT 

Figure 7 shows the distribution of mean velocities plotted 
in the form of “friction velocity” u/Ur against “friction- 
distance parameter” yU T /v . In this form the profiles are 
independent of the Reynolds number U 0 dlv and follow Von 
Karman’s logarithmic law: 

A log 

Ur V 

The constants A and B are 6.9 and 5.5, respectively. 
Comparing the present measurements with previous ones it 
is to be noted that because of the very smooth wall surfaces 
used, the value B is larger than that of other investigators. 
Making allowances for this effect, measurements are in 
fairly good agreement with those of Donch (reference 4) 
and Reichardt (reference 9). Nikuradse’s channel results 
(reference 5) differ from all other existing experiments. It is 
pertinent to mention that for the Reynolds number range in 
question (/?< 100,000) the slope of the logarithmic mean- 
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velocity distribution (the constant A) both in pipes and 
channels is larger than that established by Nikuradse’s pipe 
measurements (A = 5.75) . 

The Reynolds number has a definite influence on the 
velocity fluctuations also. Along most of the cross section 
where the influence of the viscosity is negligible the fluctua- 
tions u' , 7', and w' decrease slowly with increasing Reynolds 
number. 

Very close to the wall, according to Prandtl’s hypothesis, all 
velocities of the form Velocity/ U T must be a function of the 
friction-distance parameter only. It has already been 
pointed out that the mean velocity obeys this similarity 
law. Figure 26 shows the distribution of u' /U T as a function 
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Figure 26.— Logarithmic representation of velocity fluctuations u’. 


of yUr/v for the different flow conditions investigated. The 
following remarks can be made with reference to this figure: 

(1) For values of yU T /v<^lOO } values of u' /U T correspond- 
ing to the various cases indicate similar behavior. They 
reach a maximum at about yUr/v ~\7 . 

(2) It is believed that the similarity is actually more 
complete than indicated in figure 26. Because the micro- 
scale is not known in regions very close to the wall no length 
corrections could be applied. These corrections would of 
course be appreciably higher for the two high-velocity flows 
(jK= 61 ,600, 2d=5 in. and Ii= 12,200, 2d= \ in.) and would 
therefore bring the various distributions of u' /U T closer to- 
gether in the region in question. 

(3) The effect of viscosity is more pronounced on the 
fluctuating quantities than on the mean velocity. 

(4) Taylor pointed out in 1932 (supporting his arguments 
by Fage and Townend’s ultramicroscope measurements) that 
u' /u and w' /u approach a finite value at the wall (reference 
18). It follows from the similarity law that this value 
should be an absolute constant independent of the Reynolds 
number. Figure 10 indicates this to be true, the constant 
being (S 7 ^)*,=o~ 0 . 18 . 

It is of considerable interest to discuss the variation of the 
scale and microscale with Reynolds number. For flows be- 
hind grids where the turbulence is isotropic the scale is in- 
dependent of the mean velocity and depends on the mesh 
size of the grid. Similar behavior was found for the channel 
flow. Figure 21 shows the distributions of L y and L z for 
different velocities. These distributions indicate no con- 


sistent variation with velocity. F urthermore, measurements 
in a 1-inch channel give a value for L z five times lower and a 
ratio for L y somewhat larger than the values obtained in 
the present investigation. 

The variation of X depends, of course, on the velocity and 
channel width. The values of X decrease with increasing 
velocity; however, the variation of X with channel width is 
less than that of L. 


FULLY DEVELOPED CHARACTER OF TURBULENCE 


The flow in the channel is called fully developed if the 
variations of the mean values of the velocity and the mean 
squares of the velocity fluctuations with x are very small. 
That the mean velocity profile does not vary downstream is 
evident from the pressure-gradient measurements (fig. 19). 
The gradient in x of u' 2 was measured on the axis of the 


channel. It was found that u' 2 was indeed decreasing with 


The gradient, however, was very 
1 dp 
p dx' 


x. 


witn - ~ 


small 


as compared 


du' 1 

dx 


0.01 


I 

p dx 


Hence for all practical purposes du' 2 /dx can be neglected. 
No measurements have been made concerning du'v'/dx, 
since the scatter in the values would cover any effect. 
However, there is little doubt that du'v'/dx is of the same 
order as du' 2 /dx and that the use of equations (3a) and (3b) 
is therefore justified here. 


ENERGY BALANCE IN FLUCTUATING FIELD 


The energy equation for a two-dimensional channel has 
the form given by equation (8) and is valid throughout the 
cross section of the channel with the exception of a small 

( 

region near the wall. The term r on the left side of equa- 
tion (8) corresponds to the energy produced by the shearing 
stresses and it can be obtained directly from the measure- 
ments of r and from the mean velocity profile. The second 


term on 


the right expresses the amount of 

energy that is being dissipated because of the breaking 
down of the larger eddies to smaller ones. The term may 
be written explicitly 




T 


© ,+ (TO' + (TO 


dr'X 2 i /dw 


dw' 

,*y 


dw '\ 2 


The problem is to express these functions in terms of 
easily measurable quantities. In the case of isotropic 
turbulence Taylor solved the problem by introducing the 
microscale of turbulence X and obtained for the dissipation 


W= 1 5p 


u” 

X 2 


( 10 ) 
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It is attempted here to carry over his analysis to the ease 
of channel flow. The following assumptions have to be 
made in this connection: 

(a) The gradient of the velocity fluctuations is small in 
all three directions. With the exception of the region near 
the wall this is justifiable from the measurements. 

(b) Since only correlations of u' have been measured, it 
is assumed that the correlations of w' as functions of Z 
(Von Karman’s /-function) have the same curvature at the 
origin (Y~7= 0) as the correlation of u' has as a function 
of X) that is, 

(b 2 R z w \ 1_ 

\ dZZ ) o"A dX 2 / 0 X . 2 

Furthermore the correlations of v f as functions of X and Z 
and the correlations of w' as functions of A" (Von Karman’s 
(/-functions) have the same curvature at the origin as the 
correlation of u' has as a function of Z; that is, 



Finally, 

V dF 2 Jo \ £>F 2 Jo v C )Y 2 Jo X, 2 
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With these assumptions I lie derivatives of the fluctuations 
could be expressed in terms of the measured values of X r , 
X 77 , and X 2 : 



2 u' 2 

^ K 2 


2 u' 2 2v' 2 2v' 2 2v' 2 

l^ + T z r+ lQ r+ \ z 2 



At the middle of the channel W turned out to 
per cubic centimeter per second for //=30,800. 
isotropic relation 


2 w' 2 w' 2 l 

f V + X* 2 J 

be 2.62 ergs 
Using the 


TJ7 ,r U i z U U'~ 

W= 15/X T~2=15fJL -~2=15 m T“2 

a y Ax A 2 

the values 2.88, 4.38, and 2.24 ergs per cubic centimeter per 
second were obtained depending upon whether X r , X ?/ , or X 2 
was used. (For the value of u' 2 the algebraic mean of the 
squares of the fluctuations was used.) 

Figure 27 shows the distribution of W in the center region 



V/d 

Figure 27— Turbulent-energy balance in center region of channel. 72=30,800. 


of the channel. Taylor obtained a similar distribution of 
the dissipated energy across the channel (reference 20); 
however, his numerical magnitudes are too high since he 
substituted in the isotropic relation, in equation (10), the 
values of X y , which turn out to be less than those of X r 
and X 2 . 

From the known distributions of the energy production 
and dissipation the diffusion of energy is easily calculated 
from equation (8). It should be pointed out that because 
of the approximations involved in the calculation of the 
dissipation and diffusion terms, the conclusions derived 
from them are more or less of a qualitative nature. It is 
seen from figure 27 that at y/d~0 . 7 the diffusion term is 
zero. It was also pointed out earlier that in this region 
the /^-correlations show a considerable change in shape 
indicating a shift in energy from the lower to higher 
frequencies of the velocity fluctuations. These two facts 
suggest that the energy diffusion is associated mainly with 
the low frequencies of the fluctuations. 

The equation expressing the balance of the three forms of 
energy furnishes the following picture of the turbulent flow 
field in the channel where viscous dissipation is still neg- 
ligible: Two planes passing through points where the dif- 
fusion of energy vanishes divide the channel flow into three 
parts. The middle region, the width of which is of the order 
of Lx, receives most of its energy by diffusive action and this 
energy is dissipated here at a constant rate. In the two out- 
side regions all three energy terms increase rapidly, the pro- 
duction term being the dominant one, and their interaction 
is more involved. 

This picture of the flow field is only of a descriptive nature. 
The purpose of further investigations should be to obtain 
information on the development and mechanism of such 
an energy balance. 

LOCALLY ISOTROPIC CHARACTER OF TURBULENT CHANNEL FLOW 

The concept of locally isotropic turbulence obtained by 
Kolmogoroff requires the smaller eddies in turbulent flow 
to approach isotropy. Smaller eddies are the ones with 
length dimension l small compared with the scale of turbu- 
lence L. The smallest characteristic length in turbulent 
motion is Kolmogoroff’s rj defined as 


where e is the total dissipated energy. Clearly, to approach 
locally isotropic conditions it is necessary that 

L»77 

It follows from this hypothesis that beyond sufficiently 
high frequencies (of the order of u' /rj, say) no correlations 
exist between the components of the velocity fluctuations. 
Figure 28 shows the measured spectrum of u' 2 as compared 
with the spectrum of u'v' at y/d=QA. It is seen that 
the shear spectrum tends to zero at a frequency of about 
1500 cycles per second while F^in) still has a detectable 
value at 5000 cycles per second. This result verifies Kol- 
mogoroff ’s assumption. It should be mentioned that the 
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Fictrk 28. — Comparison of spectra of u'- and «V at y/d=0A. 


absolute values of F~^(n) are not so accurate as those of 
F^(ri) since they represent the small difference between 
two large values of hot-wire signal. 

It is to be noted that even though local isotropy was 
shown to exist in the channel flow, the values of the various 
vorticity terms (u 2 /\ 2 ) differ appreciably, indicating that the 
relative magnitudes of these terms do not constitute a 
sensitive test for the existence of local isotropy. This is 
evident from the comparison of the distributions n 2 F^{n) 
and ri 2 F~{n) (fig. 28). It is apparent from the figures 

/» 00 

that a large part of the contribution to n 2 FZJ*(n) dn 

Jo 


comes from a frequency well below 1500 cycles per second. 
Local isotropy, according to the shear spectrum, exists only 
above a frequency of 1500 cycles per second. 

CONCLUDING REMARKS 

The measurements presented here confirm the general 
conceptions concerning the mean velocity profile in a turbu- 
lent channel. The extent of the laminar sublayer, the 
velocity profile in the sublayer, and the over-all velocity 


distribution as measured here are in good agreement with 
general theoretical expectations. 

Measurements of the turbulent field show that the hot-wire 
technique is well enough developed to give consistent results 
for the intensities and correlation functions. 

Detailed measurements of the velocity fluctuations in t ho 
direction of the flow u f were carried out well within the 
laminar sublayer. It was found that the similarity law for 
u'/U r) where U T is friction velocity, holds fairly well in the 
vicinity of the wall and as a consequence the magnitude of 
u' /u, where u is local mean velocity, approaches an absolute 
constant at the wall (approx. 0.18), the value being inde- 
pendent of the Reynolds number. The magnitudes of the 
velocity fluctuations normal to the flow in the lateral and 
vertical directions v' and w' are nearly the same in the 
middle region of the channel, w' increasing more rapidly 
toward the wall. 

The measured microscales of turbulence X x , \ !n and X 2 con- 
sistently show a maximum at y/d~ 0.7, where y is lateral 
distance and d is half width of the channel; they increase 
proportionally with u ' , indicating a constant rate of energy 
dissipation W in the center portion of the channel since 
ILoc ^/X 2 . 

The scales of turbulence L y and L z in the center region are 
independent of Reynolds number and depend only on the 
channel width. The microscales, however, show a depend- 
ency on the Reynolds number. 

From the calculated magnitudes of the energy produced 
by the shearing stresses and of the dissipated energy a 
descriptive picture of the energy diffusion in the center region 
of the channel is obtained. 

The spectrum measurements of the ^-fluctuations tend 
to indicate that F^{n) behaves as n ~ 7 over a large band in 
the high-frequency region (where F^(n) is fraction of 
turbulent energy associated with band width dn and n is 
frequency). From the comparison of the spectra of u ' 2 and 
of the turbulent shear the existence of local isotropy in the 
channel flow is verified. 

Guggenheim Aeronautical Laboratory, 

California Institute of Technology, 

Pasadena, Calif., September 1 , 1949. 
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